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Abstract 

Flow and transport are studied in transparent model fractures with rough comple- 
mentary self-affine walls with a relative shear displacement u. The aperture field 
is shown to display long range correlations perpendicular to u: for flow in that 
direction, the width and geometry of the front of a dyed shear-thinning polymer 
solution displacing a transparent one have been studied as a function of the fluid 
rheology and flow rate. The front width increases linearly with distance indicating 
a convection of the fluids with a low transverse mixing between the flow paths. The 
width also increases with the flow-rate as the fluid rheology shifts from Newtonian 
at low shear rates 7 towards a shear thinning behaviour at higher 7 values. The 
width also increases with the polymer concentration at high flow-rates. These re- 
sults demonstrate the enhancement of the flow velocity contrasts between different 
flow channels for shear thinning fluids. The relative widths at low and high 7 val- 
ues for different polymer concentrations are well predicted by an analytical model 
considering the fracture as a set of parallel ducts of constant hydraulic apertures. 
The overall geometry of the experimental front geometry is also predicted by the 
theoretical model from the aperture map. 
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1 Introduction 



Transport and flow in porous media and fractured rocks are encountered in 
many engineering fields [1] and complex fluids such as polymer gels or surfac- 
tants are often involved. Applications include enhanced oil recovery (EOR), 
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drilling muds and heavy oil recovery. In EOR, for instance, polymer flood- 
ing reduce viscosity driven instabilities (a polymer solution is injected in the 
reservoir and followed by a water flood). When these complex fluids have shear 
thinning properties, experimental flow measurements (see [2]) display specific 
features such as permeability enhancements (or reductions of the effective vis- 
cosity) compared to the case of Newtonian fluids. 

These effects may be strongly influenced by fractures which are frequently 
encountered in many reservoirs and generally display a broad range of charac- 
teristic length scales. While it is customary to visualize the fractures as parallel 
plates separated by a constant distance [3J, this representation is rarely accu- 
rate: fracture wall surfaces are indeed rough and do not perfectly match [I]. 
This creates voids of various size resulting in spatial heterogeneities of the flow 
field [5116] . 

The objective of the present work is to analyze experimentally and analytically 
these velocity contrasts for shear thinning solutions flowing in transparent 
models of single fractures with rough walls and to determine their dependence 
on the fluid rheology and on the flow velocity. The experiments have been re- 
alized in a configuration in which flow is strongly channelized as is frequently 
the case in actual fractures [jQ: this will be shown to allow for analytic pre- 
dictions of the relation between the flow distribution and the apertures and, 
also, of their dependence on the rheological characteristics of the fluids used 
in the experiments. 

We have sought particularly in this work to reproduce the roughness of natural 
fractured rocks which is characterized by a broad distribution of the charac- 
teristic length scales [H] . More precisely, these surfaces can often be considered 
as self-affine [10], this means that they remain statistically invariant under the 
scaling transformation: 

h(\x,\y) = \<h(x,y), (1) 

where h(x, y) is the surface height and ( is the roughness or self-affine ex- 
ponent. For most materials including granite, £ is close to 0.8 [IT] but it is 
close to 0.5 for materials such as sandstone and sintered glass beads [T2~][T3] . 
Many experiments suggest that ( is independent on the orientation of profiles 
measured on the surface with respect to the direction of crack propagation (a 
slight anisotropy has however been recently observed experimentally on some 
materials |14j). 

The rough surfaces used in the present work are transparent milled plexiglas 
plates with an isotropic self-affine geometry of characteristic exponent ( = 0.8: 
they allow for optical flow observations by means of dyed fluids (practically, 
a transparent solution is displaced by a dyed one and the geometry of the 
front is determined by image analysis). For each fracture, two such comple- 
mentary surfaces are realized and match perfectly when brought in contact: in 
the model, both a spacing normal to the mean fracture surface and a relative 
lateral shift u are introduced in order to create a mismatch and to obtain a 
variable aperture field [T5] . 
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While the surfaces are isotropic, previous laboratory measurements and nu- 
merical investigations [71ll6fl7fl8lfT9] show that the lateral shift introduces an 
anisotropy of the permeability which is highest in the direction perpendicular 
to u. More precisely, flow channels perpendicular to u and with a length sim- 
ilar to the model appear as shown in a previous work [8]. As a result, for flow 
perpendicular to u, the overall geometry of the displacement front of a fluid 
by another of same rheological properties is well reproduced by modelling the 
fracture as a set of parallel ducts with an hydraulic aperture constant along 
the flow [S]: the present work deals exclusively with this configuration. 

The fluids used here display at low shear rates 7 a "plateau" domain in which 
they behave as Newtonian fluids of constant viscosity \i while, at higher shear 
rates, /1 decreases with 7 following a power law. Comparing the velocity con- 
trasts between the different flow paths in the two regimes allows one therefore 
to estimate the influence of the rheology since the velocity contrasts should be 
enhanced in the shear thinning case. Finally, an analytical model predicting 
the influence of this Newtonian "plateau" on the dependence of the velocity 
contrasts on 7 will be derived and compared to experimental observations. 



2 Experimental procedure 

2. 1 Characteristics of the model fracture 

The model fracture is made of two complementary rough self-affine surfaces 
without contact points: both surfaces are obtained from a transparent mate- 
rial by means of a milling machine and their size is 85 mm x 170 mm. A 
detailed description of the procedure is given in [7J: a self-affine surface h(x, y) 
is first generated numerically using the mid-point algorithm [20] with a self- 
affine exponent £ = 0.8 as observed in many materials [llj . A second surface, 
complementary from the first one, is generated and then shifted numerically 
parallel to its mean plane by 0.33 mm. The milling tool is computer controlled 
and a complex tortuous path may be imposed to obtain the self-affine geom- 
etry. Moreover, the borders of two parallel sides of the surfaces rise above the 
rough surface: they are designed so that, when they are clamped against the 
matching border of the other surface, there is a void space in the remaining 
areas. The mean planes of the facing surfaces are parallel outside these borders 
with a mean distance: a = 0.77 mm. 

The local aperture a(x,y) at a location (x,y) in the fracture plane may be 
predicted from the mathematical surface h(x, y) by the relation: 




(2) 
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where u is the lateral shift. Figure [T] shows the aperture field of the fracture 




Fig. 1. Upper figure, gray scale representation of the aperture field. Field of view: 
85 mm x 171 mm. Aperture field - mean value: a = (a(x,y))( x ^ = 0.77 mm; 

and the aperture fluctuation: a a = ((a(x,y) — a) 2 ) l ^ y ) = 0.1 mm. Shift amplitude: 
u = 0.33 mm (oriented vertically on figure). In the present work, flow is parallel 
to x direction (horizontal on the figure). Lower image: binarized aperture field with 
a threshold value equal to the mean aperture (0.77 mm). 

considered in this work: the binarized image (lower part of FigJT]) displays a 
clear anisotropy and a large correlation length perpendicular to the shift u. 
Quantitatively, this effect may be characterized by the following correlation 
function, also called semivariance [2T] : 

1 (5) = ((a(r)-a(r + 5)f), (3) 

measuring the spatial correlation of the aperture field between two points sep- 
arated by a lag vector 5. Orientations of S perpendicular (x — direction) and 
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parallel (y — direction) to the shift are of special interest. Figure [2] displays 




Fig. 2. Semivariograms normalized by 2a 2 as a function of lag distance \\5\\ (mm) 
for the aperture field displayed in FigfTJ Dotted line: correlation along the direction 
y of the shear. Solid line: correlation along the perpendicular direction x (parallel 
to the flow in the rest of the paper). 

variations of the semivariance in both directions. When the lag modulus ||<5|| 
is larger than the correlation length of the aperture field, one expects 7 to 
reach a constant value equal to 2a 2 , where a 2 = (a(x,y) — a) 2 is the variance 
of the aperture. The semivariance 7 reaches this limit, but in a very different 
way for the two orientations of 5. In the direction y parallel to the shift, 7 
becomes of the order of (and sometimes larger than) 2a 2 for ||5|| > 8 mm. 
In the perpendicular direction x, 7 never exceeds the saturation value and 
slowly increases towards it: these differences reflect the large scale anisotropic 
structure of the aperture field. Semivariograms have been computed on surface 
maps of epoxy casts of a fractured granite sample in a previous work [5j and 
display similar features: moreover, normalized curves / ~f/(2a 2 ) corresponding 
to different values of u displayed a universal variation as a function of the 
normalized lag 8/u. This suggests that results obtained in the present work 
might be extrapolated to other values of u. 

Finally it should be noted that the ratio S of the standard deviation of the 
aperture a a to the mean aperture a is only of 0.13 (Fig. [I]). This implies, 
as discussed by [22], that the fracture can be considered as " hydraulically" 
smooth with relatively small velocity contrasts between and along flow lines. 
This keeps the trajectories of the preferential flow channels relatively straight 
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and simplifies subsequent analysis. 



2.2 Experimental set-up and procedures 



The plexiglas model fracture is held vertically in a fixed position between a 
light panel and a 12 bits digital CCD camera with a high stability and dynam- 
ical range. Flow is induced by sucking a dyed solution from the top side while 
the lower side is slightly dipped into a bath containing a clear fluid. An ap- 
propriate calibration, described in reference [7J, allows one to obtain from all 
pictures of each experiment the corresponding concentration map c(x,y,t). 
Here, we focus on the geometry of the iso concentration front c/cq = 0.5 
which is determined by thresholding the concentration maps and which de- 
pends strongly on the flow heterogeneity. 



2.3 Rheological characteristics of shear-thinning soltions 




Fig. 3. Variation of the effective viscosity [i of the polymer solutions as a function 
of the shear rate 7 for two water-polymer solutions of different concentrations: 500 
ppm (□) and 1000 ppm (o). Dashed lines: Carreau functions corresponding to the 
sets of parameters of Table HJ continuous lines: truncated power law approximation. 
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In this work, we used shear thinning polymer solutions, more specifically 
water-scleroglucan solutions; they have been characterized using a Contraves 
LS30 Couette rheometer for shear rates 7 ranging from 0.016 s" 1 up to 87 
s _1 . Two different polymer concentrations equal to 500 ppm and 1000 ppm 
have been used. The rheological properties of the dyed and transparent solu- 
tions have been verified to be constant with time within experimental error 
over 3 days; the variation of the effective viscosity \i as a function of the shear 
rate 7 is displayed in Figure [31 The variation of \i with 7 is well adjusted by 
a Carreau function: 

f 1 =~ l , ^ 9 ,l-n (/4)-/Ax>)+/W (4) 

(l + (^) 2 )~ 

The values of the corresponding rheological parameters for the polymer solu- 
tions characterized in the present work are listed in Table [TJ fi^ is too low to 
be determined within the available range of shear rates (7 < 87 s _1 ) and it has 
been taken equal to the viscosity of the solvent (i.e. water). In Eq.fllJ), 70 cor- 
responds to a crossover between two behaviors. On the one hand, for 7 < 70, 
the viscosity fi tends towards the limiting value fio, and the fluid behaves as a 
Newtonian fluid. On the other hand, if 7 > 70, the viscosity follows a power law 
variation reflecting its shear thinning characteristics with fi oc ^ n ~ 1 \ For each 

Polymer Cone. n 70 /io 

ppm s" 1 mPa.s 

1000 0.26 ±0.02 0.026 ±0.004 4490 ± 342 

500 0.38 ±0.04 0.077 ±0.018 410 ± 33 

Table 1 

Rheological parameters of schleroglucan solutions used in the flow experiments. 

experiment, the flow rate is kept constant at a value between 0.01 ml/min and 
5 ml/min (corresponding mean flow velocities: 0.0003 < v < 0.14 mm.s" 1 ). 
Under such conditions the typical shear rate 7 ~ v/a ranges between 4.10~ 4 
s _1 and 0.18 s _1 . The latter value is far below the shear rate corresponding to 
the second Newtonian plateau (/i = /ioo) and this limit will not be considered 
in this work. On the contrary, the lowest values of the typical shear rate are 
much lower than 70 : the Newtonian "plateau" in the rheological curves may 
therefore have a crucial influence of the flow properties. 

In order to obtain an analytical expression accounting for the effect of the 
fluid rheology on the velocity fluctuations, the rheological law of the fluids is 
approximated in section 0] by a truncated power law. 

When 7 < 7o, the viscosity /i(7) is considered as constant and equal to /i ; 
for 7 > 7b, ^(7) is assumed to follow a power law ^(7) = m^ n ~ x in which 



7 



m = /io/7o™ 1 - The parameters n, 70 and fio are obtained from TabJT] 



3 Flow velocity dependence of front geometry 

Two important features of the displacement front have been observed: (a) its 
geometry depends on the flow rate Q, and (b) its width parallel to the flow 
increases linearly with the distance from the injection side. The first point 
is illustrated by Figure H] in which two fronts measured are compared during 
fluid displacements at two different flow rates, but for a same polymer con- 
centration (1000 ppm): the front width parallel to the flow direction is clearly 
larger at the highest flow rate. 

The broadening of the front may be characterized quantitatively from the vari- 
ation of the mean square front width, a x (t) = {(x(t) — x(t)) 2 ) 1//2 , as a function 
of the mean distance x(t) of the front from the injection side (Figure [5]). For 
all values of Q, o~ x (t) increases linearly with x(t). In the next section, this will 
be shown to result directly from the underlying channelized structure of the 
aperture field. The width of the front after a transit time t corresponds then 
directly to the product tAv where Av is the velocity difference between the 
different channels (the transverse exchange between channels is too small to 
allow one to reach a diffusive spreading regime). 

At all distances, the width a x (t) increases with the flow rate Q but with a 
particularly sharp variation between Q = 0.1 ml/min and Q = 0.5 ml/min. 
It will be seen that, at this transition flow-rate, the shear rate at the fracture 
walls becomes of the order of 70 (the threshold value above which the fluids 
display shear thinning characteristics). 



4 Modelization 

4-1 Flow between parallel plates with a constant gap 

We compute the flow of the polymer solutions between parallel plates using 
the same approach as in reference [23]. The relation between the longitudinal 
pressure drop and the velocity profile in the gap is obtained by using the trun- 
cated power law model discussed in section [2731 The flow field is unidirectional 
and parallel to x so that v x (z) is the only non zero velocity component. The 
strain rate is given by: j(z) = dv x /dz. 

At low pressure gradients, the fluids behave like a Newtonian fluid with a 
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Fig. 4. Displacement fronts at two different flow rates, (a) Q = 1 ml/min 
(v/v c ~ 100); (b) Q = 0.1 ml/min (v/v c ~ 10). Polymer concentration: 1000 ppm. 
Vertical mean flow. Front width perpendicular to mean flow: 85 mm. Front exten- 
sion along the flow: 53 mm (top curve), 34 mm (bottom curve). Mean distance of 
the front from the inlet equal to half the fracture length. 
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Fig. 5. Variation of the mean front width a x {t) = ((x(t) - x(t)) 2 ) 1 / 2 as a function 
of the mean distance x(t) from the inlet side of the model. Polymer concentration: 
1000 ppm. Flow rates: Q = (A)2.0,(+) 0.5, (o) 0.1, (□) 0.02 ml/min. 

constant viscosity /i and the resulting velocity profile is parabolic and sym- 
metrical between the walls. Then, the shear rate is zero half way between the 
fracture walls and reaches a maximum at their surface where 7 = 6v/a (a is 
the distance between the plates). This value of 7 is proportional to the mean 
flow velocity or, equivalently, to the pressure drop. 

As the flow rate keeps increasing, 7 becomes larger than 70 and the non New- 
tonian characteristics of the fluid modify the velocity profile. The mean flow 
velocity v c corresponding to the transition between the two regimes satisfies: 
v c = 070/6 and the corresponding pressure gradient is: VP C = 2/io7o/ a - 
As v increases above v c , the layer where the shear rate is higher than 70 be- 
comes thicker and the velocity profile v x (z) is no longer parabolic: the full ex- 
pression may be derived analytically and is given in Eq. (5) of reference |23j . 
The mean velocity, v, can then be computed by integrating v x (z) over the 
fracture gap, leading to: 



" 2 -VP. 



12(2ri + l)/i 



VP,_ 3 „ ,VP 



n){— — )~ 6 + 3n 



1-n 



VP/ v vp 



(5) 



We consider now the case of shear thinning fluids such that n > and 
(1 - n)/n > -1. Then, the leading term in Eq.© is (VP/VP c ) (1_n)/n and, 



10 



therefore, when VP >> VP C , Eq. (jSj) becomes: 

a 2 fVP\" 

where /i e // = fio {^io/ a ) 1 ~ n ((2n + l)/3n) n . This is similar to the generalized 
version of Darcy's law often applied to the flow of non Newtonian and, more 
specifically, to power law fluids in porous media [24|26f27j . 



4-2 Flow in rough fractures 



In this part, we focus on the variations of the velocity in the plane (x, y) of the 
fracture and we assume therefore a two-dimensional flow field v(x, y) equal to 
the average of the fluid velocity profile over the gap with v(x, y) = (v(x, y, z)) z . 
The development with time of the front (represented by the iso concentration 
lines c/cq = 0.5) will now be analyzed by assuming that its points move at the 
local flow velocity v(x, y) and an analytical model predicting the global front 
width will be developed. 

This model is based on the results of a previous work [5] demonstrating that, 
in such systems, the aperture field is structured into channels perpendicular to 
the lateral shift u of the surfaces. For a mean flow parallel to these channels, the 
paths of the tracer particles have a weak tortuosity; also, the velocity variations 
along these paths are small compared to the velocity contrasts between the 
different channels. Under these assumptions, the velocity of a particle located 
at a distance y, perpendicular to the mean velocity, satisfies: 

v(x,y) » v(y)e" x , (7) 

where e x is the unit vector parallel to the mean flow. Note also that, in the 
geometry discussed in this section, there are no contact points between the 
walls of the fractures: this avoids to take into account the large tortuosity of 
the flow lines in their vicinity. 

If the fluid is Newtonian with a constant viscosity, then, for each channel, the 
velocity is related to the pressure gradient VP by relation (jB]) with n = 1; 
a is now an equivalent (or hydraulic) aperture associated to each channel 
and noted a(y) and the equation represents the classical linear equivalent of 
Darcy's law for fractures. Previous studies have shown that, for relatively small 
aperture fluctuations, this hydraulic aperture is well approximated by the ge- 
ometrical aperture [5|2"2"] : this suggests that a(y) 2 can be taken equal to the 
mean of the average of the square of the local apertures along the direction 
x i.e. a(y) 2 = (a(x,y) 2 ) x . The validity of this assumption has been tested 
numerically previously for a similar geometry [8] in the case of a Newtonian 
fluid: these simulations used the lattice Boltzmann method to determine the 
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2D front geometry at all times: except for fine scale details, the profile x(y, t) 
of the distance of the front from the inlet at a given time t follows very closely 
the variations of a(y) 2 . 

For a power law fluid such that n < 1, the velocity satisfies the non linear 
generalized relation (Q. We seek now to generalize to this case the relation 
between the front geometry and the aperture variation established for the 
Newtonian fluids: the aperture field is still assumed to be strongly correlated 
in the flow direction, allowing one to consider the fracture as a set of parallel 
ducts. 

We consider particles starting at t — from the inlet of the model at differ- 
ent transverse distances y and moving at different velocities v(y). Then the 
distance x of the particles from the inlet at time t after the injection satisfies 
x(y,t) = v(y)t so that the mean distance of the front from the inlet side is 
x{t) = (v(y)) y t = vt and: 

Moreover, the mean square deviation a x (t) = ((x(y,t) — x) 2 ) 1 / 2 should satisfy: 
&x{t) = <y v t where a v is the mean square deviation of the velocities in individual 
channels from their mean value v. Combining the previous relations leads to: 

x(t) v ' 1 ' 

This equation shows that there is a direct relation between the front geome- 
try and the variations of the velocity from one channel to another: for power 
law fluids, the latter are related to the variations of the hydraulic aperture by 
Eq. (jSJ). In order to estimate these variations, we introduce a modified reduced 
aperture deviation Sh defined as the ratio between the standard deviation of 
the hydraulic aperture a(y) to its mean. The parameter Sh is equivalent to 
the reduced aperture deviation S defined in section 12.11 but the geometrical 
aperture is replaced by the hydraulic one. Here, we are interested in weakly 
fluctuating systems, i.e. for which both S and Sh are small compared to one. 
In addition, the hydraulic aperture a(y) is observed to follow a Gaussian dis- 
tribution. Moreover, Eq. (jSJ) shows that, for a given pressure gradient VP, 
v scales as a n+1 / n : together with the above assumptions, this leads to the 
following relation between the reduced velocity fluctuations a v /v and Sh'- 

^ = ^S h . (10) 

v n 

Combining Eqs. (Q and (jTUI) . leads to: 

— 77T- = Oh- I 11 ) 

x(t) n 
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5 Quantitative comparison between the experiments and the model 

In the present experiments, the polymer solutions are expected to behave like 
Newtonian fluids as long as the shear rate 7 is everywhere lower than the 
critical value 70 (see Table [T|). As the flow rate increases, the critical shear 
rate 70 is first reached at the wall of the fracture where 7 is highest. If the 
fracture is modeled as two parallel plates separated by the mean aperture a, 
then 7 = 70 at the walls when the mean flow velocity is v c = cl^q/Q. Above 
this velocity, the shear thinning properties of the fluids influence the flow and 
enhance the velocity fluctuations. 




Fig. 6. Experimental variation of the ratio a v /v = a x /x in the model fracture 
(Sh ~ 0.02) as a function of the normalized mean velocity v/v c for shear thinning 
solutions of concentrations 1000 ppm (o) and 500 ppm (□). Horizontal dotted lines: 
theoretical values of a v jv computed from Eq. flTl) for a Newtonian fluid (n = 1) 
and for 500 ppm (resp. 1000 ppm) polymer solutions (n = 0.38, resp. 0.26). Solid 
and dashed lines: variations of a v /v a as function of v/v c computed by integrating 
Eq. §J) 

Fig. [6] displays the experimental variations of the normalized velocity fluctu- 
ations (equal to the normalized front width a x (t)/x(t)) as a function of the 
reduced velocity v/v c for both polymer solutions. The values predicted by 
Eq. ( TTTT) for a Newtonian fluid (n = 1) and for power law fluids with the same 
index as the two solutions are also plotted. 

For v/v c < 1 the experimental values are similar for both solutions and close 
to the theoretical prediction for n = 1 (horizontal dashed line). For v/v c ^> 1, 
a x {t)/x tends toward values of the order of those predicted by Eq. fTTTj) and 
increasing with the polymer concentration. 
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Eq. (fTTj) provides therefore a good estimate of the velocity fluctuation inside 
the fracture both for low, i.e. v/v c < 1, and high flow rates corresponding 
to v/v c >> 1. The increase of a x (t)/x between the Newtonian and shear 
thinning regimes and also, at high velocities, with the polymer concentration 
confirms the enhancement of the velocity contrasts between the channels for 
shear-thinning fluids. 

Between the limiting values v/v c < 1 and v/v c ^> 1, fluid velocity variations 
within the fracture may be estimated by applying Eq.([5]) in each channel (as- 
sumed to be of constant hydraulic aperture): this equation takes into account 
the coexistence in the fracture gap of layers where the fluid has Newtonian and 
non Newtonian properties. The normalized velocity fluctuations a v /v obtained 
by these computation are displayed in Figj6] for the two polymer concentra- 
tions together with the experimental variations of the normalized front width 
a x (t)/x(t). 

In agreement with the theoretical curves, a x (t)/x(t) starts to increase when 
the velocity v becomes larger than v c (v/v c > 1) for both polymer solutions. 
However, although the limiting value for v/v c ^> 1 is the same as predicted, 
the increase of a x (t)/x(t) above v c is slower than expected: actually, the the- 
oretical predictions represent an upper bound for the observations. 
This difference may be due in part to the use of a simplified version of the 
rheological curve displaying a transition sharper than the actual one between 
the Newtonian and shear thinning regimes. Also, the aperture of the parallel 
channels introduced in the model is assumed to be constant: this also leads 
to a transition between the Newtonian and power law regimes which is faster 
than the actual one. 

A step further in the interpretation is the comparison of the experimental 
shape of the fronts with that estimated from the channel model. In section |4T2| 
the normalized distance x(y,t)/x(t) of the front from the inlet has been pre- 
dicted to be equal to v(y)/v (see Eq. [S]). An experimental front profile normal- 
ized in this way is plotted in Fig. [7] as a function of the transverse distance y 
together with the variation of the theoretical normalized velocity v(y)/v. The 
velocity v(y) is estimated from Eq. [6] in which the aperture a is replaced by 
the mean value a(y) defined in section [4721 

The most remarkable observation is the fact that both the experimental and 
theoretical fronts have not only the same width but also nearly the same ge- 
ometry. These results are very similar to those of numerical simulations for 
Newtonian fluids [8]: they demonstrate the validity of the generalization in 
Eq. ([6]) of the Newtonian model. Fine scale details predicted by the theoreti- 
cal model are however not observed in the experimental front: this difference 
may be due in part to viscous drag forces between parallel layers of fluid mov- 
ing at different velocities in the fracture plane. These forces may smoothen the 
local velocity gradients and rub out small scale features of the front without 
changing the large scale velocity variations: this results in a bumpy front with 
a typical width of the structures of the order of 10 mm. This latter value is 
of the order of the correlation length in the direction perpendicular to the 
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Fig. 7. Thick solid line: Experimental normalized front profile x(y, t)/x(t) as a func- 
tion of the transverse distance y (mm) for v/v c = 200 for a 1000 ppm shear thinning 
polymer solution. Dotted line: theoretical variation of the normalized velocity v(y) jv 
in the parallel flow channel model. The front is displayed just before the displacing 
fluid starts to flow out of the fracture. Dashed lines: characteristic deviations of the 
distance x(y) from its mean value x. 



6 Discussion and conclusions 

In the present work the enhancement of velocity fluctuations for shear thin- 
ning fluids has been studied in a single fracture with rough, self-affine walls. 
The two wall surfaces are perfectly matched and are positioned with both a 
normal and a lateral shift. This results in an anisotropic aperture field well 
characterized quantitatively by the semivariograms of the aperture both in 
the direction of the shift and perpendicular to it. The characteristics of these 
semivariograms are in agreement with previous experimental measurements 
on granite samples [SJ. Parallel to the shift, the aperture field is correlated 
over a distance of the order of 10 mm above which the value of the normal- 
ized semivariograms is of the order of 1. In the other direction, the correlation 
subsists over the full fracture length. 

This observation has allowed us to model the fracture as a set of parallel ducts 
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perpendicular to the shift and with an hydraulic aperture constant along their 
length. These assumptions lead to specific predictions on the dependence of 
the width and of the geometry of the front on the velocity of the fluid and its 
rheology: these predictions deal with the case of a mean flow parallel to the 
channels which was the configuration used in the present experiments. This 
model generalizes a previous one developed for Newtonian fluids and which 
has been validated by numerical simulations [8]. 

The variation of the front width with the velocity could first be predicted. At 
low flow rates, the viscosity of the solutions is constant (Newtonian "plateau" ) 
but non Newtonian effects become important for faster flows: this results in 
an increase of the velocity fluctuations -and of the front width. This variation 
occurs when the shear rate at the fracture wall becomes larger 70, i.e. the 
shear rate corresponding to the crossover between the Newtonian plateau and 
the power law regimes: 70 is reached for a mean flow velocity v c = 070/6. At 
still higher flow velocities of the order of 100 x v c , both the normalized veloc- 
ity fluctuations and the normalized front width reach a new constant value 
with a good agreement between the experimental results and the theoretical 
expectations. 

These results validate the prediction of an enhancement of velocity contrasts 
for shear thinning channelized flows in fractures. The experimental increase of 
the front width with the mean velocity v right above the threshold value v c is 
however slower than the predictions. The origin of this discrepancy might be 
investigated by a more refined theory taking into account the full rheological 
characteristics of the fluid and aperture variations along the flow. 
The theoretical model also allows to reproduce well the experimental front 
geometry for length scales larger than 10 mm which represents the correlation 
length of the aperture field in the direction perpendicular to the flow channels. 
Future work should investigate the influence of transverse velocity gradients 
on the shape of the front for different types of fluids. 

The results obtained in the present work demonstrate therefore clearly that 
approaches developed to analyze channelized Newtonian flows in fractures can 
be generalized to non Newtonian fluids and allow to predict, for instance, the 
variation of the velocity contrasts with the rheology. 

Numerical studies in 2D networks [26.27J had similarly shown that the flow of 
shear thinning fluids is localized in a smaller number of preferential paths than 
for Newtonian ones. It has been suggested that these effects might account for 
the permeability enhancement for such fluids mentioned in the introduction: 
the results obtained in the present paper may therefore be usefully applicable 
to the numerical simulation of non Newtonian flows in fracture networks. 
A number of questions remain however open and need to be considered in 
future studies. First, the present experiments have been realized with a mean 
flow parallel to the channels created by the relative shift of the wall surfaces. 
It will be important to compare these results with the case of flow perpen- 
dicular to these channels: velocity fluctuations in the directions parallel and 
perpendicular to the flow should then be significantly different from those in 
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the present experiments. Eq. [5] should, for instance, be modified. Also, the 
present experiments deal with relatively short path lengths such that trans- 
verse exchange between channels may be considered as negligible: the results 
obtained may therefore be different for longer path lengths. It is also possible 
that the spatial correlation of the velocity field will eventually decay at very 
long distances although this has not been observed in our experimental model. 
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